Two-Gaussian excitations model for the glass transition 
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We develop a modified "two-state" model with Gaussian widths for the site energies of both 
ground and excited states, consistent with expectations for a disordered system. The thermodynamic 
properties of the system are analyzed in configuration space and found to bridge the gap between 
simple two state models ( "logarithmic" model in configuration space) and the random energy model 
("Gaussian" model in configuration space). The Kauzmann singularity given by the random energy 
model remains for very fragile liquids but is suppressed or eliminated for stronger liquids. The sharp 
form of constant volume heat capacity found by recent simulations for binary mixed Lennard Jones 
and soft sphere systems is reproduced by the model, as is the excess entropy and heat capacity of a 
variety of laboratory systems, strong and fragile. The ideal glass in all cases has a narrow Gaussian, 
almost invariant among molecular and atomic glassformers, while the excited state Gaussian depends 
on the system and its width plays a role in the thermodynamic fragility. The model predicts the 
existence of first-order phase transition for fragile liquids. The analysis of laboratory data for 
toluene and o-terphenyl indicates that fragile liquids resolve the Kauzmann paradox by a first-order 
transition from supercooled liquid to ideal glass state at a temperature between T g and Kauzmann 
temperature extrapolated from experimental data. We stress the importance of the temperature 
dependence of the energy landscape, predicted by the fluctuation-dissipation theorem, in analyzing 
the liquid thermodynamics. 



I. INTRODUCTION 



In the search for understanding of the glass transition 
phenomenon, attention has been focused overwhelming 
on the dynamic aspects of the behavior of supercooling 
liquids, 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 This is natural in view of 
the general agreement that it is the falling out of equilib- 
rium, at a temperature that depends on the cooling rate, 
which provokes the observed "drop" in heat capacity at 
T g . In other words the glass transition phenomenon ob- 
served experimentally is an entirely kinetic phenomenon. 
However, this approach leaves unresolved a basic ther- 
modynamic question that has troubled glass scientists 
for the best part of a century. 

The thermodynamic problem concerns the course of 
the entropy in excess of that of the crystal (or any other 
state whose entropy vanishes at K) during cooling of the 
equilibrated liquid state. First posed in 1930 by Simon 
for the particular case of glycerol^ and after for a vari- 
ety of substances by Kauzmann^ the question concerns 
what physical process occurs to avoid the liquid entropy 
intersecting that of the crystal, as simple extrapolation 
of the observed entropy changes with decreasing tem- 
perature would require for all fragile liquids^ Unless it 
can be shown generally that the liquid becomes mechan- 
ically unstable during cooling, (hence has no option but 
to crystallize), the resolution of this problem requires a 
thermodynamic description of the liquid entropy which is 
independent of equilibration time scales. A mechanical 
instability due to the vanishing of the nucleation bar- 
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rier was Kauzmann's resolution^ of what has become 
known as the Kauzmann paradox (kinetic phenomenon, 
T g , avoiding a thermodynamic crisis, at Tk). Although 
this resolution has been given recent support from cer- 
tain crystallizable spin-glass model studies^ there is a 
broad belief that the Kauzmann paradox demands a more 
general resolution. 

While there have been a number of insightful in- 
vestigations of the thermodynamic properties of glass- 
formers, using the configuration space energy landscape 
approachji£ii2i2& there have been surprisingly few at- 
tempts to provide theoretical functions to describe the 
liquid thermodynamics in terms of underlying models. 
Early attempts were focused on polymers for which quasi- 
lattice models were plausible. Considering the case of at- 
actic polymers, for which no low energy crystalline state 
exists, Gibbs and Dimarzk»2I argued that a thermody- 
namic (equilibrium) transition of second order, at which 
the configurational entropy vanishes, must set the limit 
to supercooling of the liquid state of the polymer. It has 
been broadly supposed that a similar transition might 
apply to liquids^ 2 , though, without the polymer basis, 
there is so far less theoretical justification for this. Fur- 
thermore, Stillinger— has argued that such a transition 
is not possible in principle, though how closely such a 
transition could be approached has not been discussed. 
By contrast, a free volume model by Cohen and Grest^ 
has suggested that, ideally, the transition to the ground 
state glass should be of first order, though no experimen- 
tal example has been identified. On the other hand, spin 
models^ and their application to coarse-grained models 
of dynamics in structural glasses f 2 4 *2 5 . have treated the 
thermodynamic component of the problem as trivial, to 
be resolved by the thermodynamics of uncorrelated exci- 
tations. 
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Indeed, it has been long known that simple un- 
corrected excitation (or defect) models of amorphous 
solids, 26 i 27 i 28 i 29 can give a good account of the entropy- 
temperature relationS 9 *^ particularly in elemental cases 
like selenium^ 9 " These show that the Kauzmann limit 
paradox can result as a consequence of an unjustified ex- 
trapolation of the entropy vs temperature relation, which 
should be continuous, though rapidly varying, in the 
vicinity of the Kauzmann temperature. Unfortunately, 
excitation models in their usual forms (in which the exci- 
tations are presumed to be non- interacting) , predict the 
occurrence of a heat capacity maximum above TgjSLS 9 . 
which in practice is only found in some strong liquids^i^ 

The simple two state model has recently reappeared 
under a new name, the "logarithmic" model based on 
its properties in configuration spaced Like its real space 
predecessors, the logarithmic model has the problem of 
predicting a heat capacity maximum where none is found 
(though when combined in configuration space with a 
Gaussian component, this problem is avoided^ see be- 
low). In a variant of such models, Tanaka— ^ has in- 
troduced a two order parameter Landau model for the 
thermodynamics of glassformers where bond length and 
orientation are distinguished. 

A defect model with behavior much like that to be de- 
scribed in this paper (despite a quite different starting 
point) is the interstitialcy model of Granato. 36,37 This 
model posits a single entropy-rich defect (the intersti- 
tial defect of crystalline metals) and obtains the temper- 
ature dependence of the defect concentration from the 
temperature dependence of the shear modulus. The co- 
operativity missing from earlier two state models, or in- 
cluded ad hocr^^ is built in through the proportion- 
ality of the defect energy to the shear modulus. The 
latter decreases strongly with temperature, leading to 
laboratory-like heat capacities and a phase transition at 
lower temperatures - which is assigned to a return to the 
crystal state. Alternatively, Wolynes and co-workers 3 ^ 9 , 
have described a mosaic model in which the inter-domain 
boundary energies play a vital role in the thermody- 
namics. With the appropriate assumptions, this model 
can resolve the Kauzmann paradox in the same way as 
does the random energy model^S^i*^ the system sim- 
ply running out of states at a singular (Kauzmann) tem- 
perature. The sudden, latent heat-free, transition to 
the ground state is described as a "random first order" 



transition* 3 ** 3 * the latent heat of the normal first or- 
der transition having been given up continuously over 
the supercooling temperature range. Both mosaic 3 *^ 
and constrained excitation* 2 * 2 *! 2 * models prove capable 
of predicting important dynamic features of glassform- 
ers, such as the decoupling of viscosity from diffusivity on 
approach to the glass transition from abovei* 3 " However, 
in this paper we are concerned only with the thermody- 
namic problem. 

Most theoretical models now gain their support from 
molecular dynamics computer simulations but, because 
of their time scale limitations, these cannot be expected 
to help much with the long time aspects of glass transi- 
tion problem. The Gaussian distributions of configura- 
tional states found in several casesai* 1 *^ 6 *^ 7 *^**?" in the 
shorter relaxation time domain (which, however, covers 
most of the inherent structure energy range between the 
extrapolated Kauzmann temperature and the high tem- 
perature limit) would imply the existence, at lower tem- 
peratures, of a Kauzmann-like singularity. Clearly some- 
thing has to change between the lowest temperature of 
these simulations and the vanishing entropy temperature, 
if Stillinger's argument is to be upheld. In the laboratory 
behavior of the closest relatives of the most simulated 
system, binary mixed Lennard- Jones (LJ) j*^* 6 **?. what 
changes is the state of the system: it crystallizes, leaving 
the problem unresolved. 

Recent simulation of the thermodynamic behavior of 
a small periodic box of the mixed soft sphere system by 
Grigera and Parisi** - Yu and Carruzzoj**^ and De Pablo 
and co-workers^ 9 - suggest, however, that if crystalliza- 
tion does not occur, and equilibrium is maintained, then 
the heat capacity continues to increase. Finally, it peaks 
sharply and decreases to zero, like a narrowly avoided 
Kauzmann singularity. Yan et al~ suppose that all pos- 
sible states of the system have been explored, though this 
is not yet proven (simulations by Yu and Carruzzoi**^ 
actually indicate that heat capacity drops due to insuffi- 
cient sampling) . In a separate study by Debenedetti and 
Stillinger j* 3 . a range of behavior intermediate between the 
simple two state model and the singularity of the ran- 
dom energy modela^ has been illustrated by adopting, 
ad hoc, an additive mixture of two-state (logarithmic 
model) and Gaussian (random energy model) distribu- 
tions. The behavior seen by Yan et al.^ and required 
by experiment** 7 " is found for Gaussian-rich mixtures. 



In the present paper we show how thermody- 
namic behavior of the sort obtained on soft sphere 
mixtures*?*^****^ can be reproduced by a modified ver- 
sion of the simple excitations model in which the sin- 
gle (or few) excitation energy(ies) of the original models 
is(are) replaced by a more physically reasonable Gaus- 
sian distribution, the centroid of which may lie near but 
generally below the value of the original excitation en- 



ergy, and the width of which may vary. The existence of 
such character has been suggested by analysis of spectral 
band-shapes for glasses and liquids*?- and recently, also, 
the vibrational density of states of glasses of different Ac- 
tive temperatures in which a quasi-two state behavior is 
found for the temperature dependence of the vibrational 
density of statesi**§* We note that the Gaussian analy- 
sis that is often used to describe the widths of spectral 
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FIG. 1: Real space (RS) energies of perfect crystal sites and defect sites (a). Real space energies and populations of the ideal 
glass and liquid sites at temperatures T\ — T g and T2 respectively (b). Configuration space (CS) location of system points at 
Ti = Tk and T2, and the excitation profile (sigmoid) (c). Quasi- Gaussian enumeration function for configurational states (d). 



bands, is only appropriate if the modes are localized - 
which is a poor approximation when dealing with the vi- 
brational density of states, even for the boson peaki£& 
Other spectroscopic evidence for distinct broken bond 
excitations in glasses has been given for weak network 
liquids 57 and recentlji^S revived in connection with the 
boson peak controversy. 



II. MODEL REPRESENTATIONS IN REAL 
SPACE AND CONFIGURATION SPACE 
(ENERGY LANDSCAPE). 

The model assumes the presence in the condensed 
phase of degrees of freedom which, in real space (RS), can 
exist in ground (low-energy) and excited (high-energy) 
statesi^ 7 . To visualize the model and its relation to crys- 
tal defect physics on the one hand, and to energy land- 
scape representations on the other, we use Fig. fl] The 
distribution of energies in the real space ground state, in 
the case of the crystal, is a delta function, on the unit 



cell length scale, and a small number of delta functions 
on the per molecule scale (Fig. QJi) . The defect states are 
likewise few in number and well defined in energy. In the 
glass, however, the sites are not all equivalent on these 
length scales, and a Gaussian distribution in energy is ex- 
pected, both for molecules in the ideal glass and for the 
elementary configurational excitation (or defect) states 
which we suppose to exist (Fig. fit). Excitations may be 
related to coordination defects in covalent materials or to 
local distortions or packing strains in molecular crystals. 

The distribution of energies ej in real space is charac- 
terized by two Gaussians Gj(ej) where i = 1 stands for 
the ground state and i — 2 stands for the excited state. 
The Gaussian function Gj(ej) is defined by the average 
eoi and the variance Oi. In addition to the change in 
energy, the creation of a local defect may result in an 
entropy increase related either to a change in the vibra- 
tional or configurational density of states^ The entropy 
change per molecule of the glass is s — ASo/Nk^. 

The potential energy in configuration space (CS) is a 
hypersurface depending on all degrees of freedom of the 
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disordered liquid. The overall configuration space is de- 
composed into basins of local potential energy minima 
termed inherent structures ^iSS The ideal glass, in config- 
uration space, is represented by the lowest energy basin 
on the energy landscape, and any excitation of defect 
states will lift the energy to one or other of the higher 
energy basins. The more defects in the real space quasi- 
lattice, the higher the energy of the configuration space 
basin occupied by the system (Fig. [TJ;) . Thus as the in- 
tensity or occupation number of the second real space 
Gaussian increases (cf. dashed to solid lines in Fig.^ 3 ), 
the system point in configuration space moves higher on 
the landscape (Fig.nj). 

The distribution of basin energies (CS) is found, by 
simulation studies, to conform to a Gaussian (Fig. QJi), 
but it is barely possible to distinguish between a Gaus- 
sian and the binomial distribution that is expected for a 
two-state system, except at the wings, which are unex- 
plored in any simulations on accessible time scales. The 
difference must diminish further for the case where the 
excitation energy is distributed, as in our model. Where 
a given excitation can occur at any energy in the real 
space distribution, the states in the configuration space 
distribution are only occupied on the low energy side of 
the maximum of the Gaussian (Fig. [TJi). The energies 
in this half Gaussian, though, are uniformly higher than 
those in the crystal manifold, which is very narrow (Fig. 
12;) , because the crystal generates very few defects before 
it becomes thermodynamically unstable and melts. 

The density of inherent structures identified with 
basins of depth (f> defines the enumeration function u)(<p). 
Following Dcrrida^ and Wolynes?Si u)(<j>) can be found 
by summation over all populations of the excited state 
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where N 2 is the number of excitations out of N molecules 
in the system, x = N2/N is the population of the excited 
state. The function C(x) in Eq. (JIJ is the number of 
realizations of a given distribution of molecules between 
the ground and excited states 
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FIG. 2: (a): a;* = x(<f>*) calculates as the root of Eq. I|21|l . (b) 
u>(4>) from Eq. 110H and its separation into the ideal mixture 
term "so" (first summand in Eq. flQ) ) and the Gaussian term 
"G" (second summand in Eq. HlUl h (c): lj(4>) from Eq. 1101 . 
The lines in (a)-(c) refer to different temperatures: T = 500 
K (solid line), T — 200 K (dashed line), and T = 10 K (dash- 
dotted line); \ 2 /k B = 200 K, Ai/fe B = 100 K, e /k B = 800 K, 
so = 2.0. The dash-dotted line in (b) shows the enumeration 
function of the two-state, "TS", model [Eq. itTSI l. 
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(2) where £q = £02 — eoi is the excitation energy and 



The distribution of basin energies P{4>, x) in configura- 
tion space [Eq. Q] can be obtained from the real space 
Gaussians (Fig.^O, 



P(<j>, x)= I 5 (0 - xe 2 - (1 - x)e 1 ) G 2 (e 2 )G 1 (e 1 )de 2 de 1 . 

(3) 

Equation @ gives a Gaussian distribution of basin en- 
ergies with the average and variance dependent on the 
population of excited states in real space 



P{4>, x) = [2-Ka] exp 



(<j) - xe ) 2 
2a(x) 2 



(4) 



a(x) 2 = {l-xf(rl + x 2 al 



(5) 



is the variance of basin energies. 

In the thermodynamic limit N — > 00, the behavior of 
the sum in Eq. is defined by its largest summand. 
One finds 



u{4>) = s(0,x(4>)), 



(6) 



where the population at maximum x((f>) is obtained from 
the stationary point 



ds((j), x) 



dx 



= 0. 



(7) 



x— x(<p) 
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The function s(<p,x) = S(4>,x)/Nk& is the entropy per 
molecule at a given population of excited states. From 
Eqs. (JTJ), ©, and |Q one finds 

s{(f>, x) = s (x) 2<t(s) 2 ' ^ ' 

where sq(x) is the entropy of an ideal binary mixture 

sq(x) — xsq — x In x — (1 — x) ln(l — x). (9) 

The requirement to maximize the entropy at a given 
basin depth makes the population x a function of (f> and, 
therefore, transforms the enumeration function 

u{4>) = so(x{4>)) — - (10) 

2cr(a;(</»)) 2 

into a generally non-Gaussian dependence on <f> (Fig. [2J . 
The solution x(<f>) is a root of Eqs. Q and © (Fig.^Ji). 

In the static energy landscape picture a system at con- 
stant volume has a unique energy landscape in the config- 
uration space that is fixed by the intermolecular potential 
for the particles in the system. From this viewpoint, it 
seems natural to assume that the variance a is indepen- 
dent of temperature. If, in addition, the dependence of 
x{4>) and a((f>) on <f) is neglected in Eq. ijTU|) . one arrives 
at the standard Gaussian model equivalent to the ran- 
dom energy model introduced by Derrida42*2i The Der- 
rida model predicts an ideal glass transition at the Kauz- 
mann temperature fceTV = a j \J 2sq(x(4>*)) (</>*) = 
at T = Tk) and a hyperbolic temperature dependence of 
the average basin energy 

4>*=a-b/T. (11) 

The dependence of the type given by Eq. (|11|) has in- 
deed been observed in several simulation studies of binary 
Lennard- Jones (LJ) mixtures^*^ although deviations 
from this law have also been reportedi^^LS^JiiiS^ The 

heat capacity per molecule cy = Cv/NttB — k^ 1 dxj>* / dT 
then varies as cv oc 1/T 2 . 
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FIG. 3: Temperature dependence of the effective excitation 
energy e(x*) [Eq. 1221 1. e(0) and e(l) denote the excitation 
energy at zero and infinite temperature, respectively. 



In the limit when the RS Gaussians are much nar- 
rower than the excitation energy, Oi eo, one gets the 
two-state modeliSi The Gaussian term in Eq. l|lUfl then 
generates a delta function in the density of states 
requiring x(<fi) = 4>/eo- This solution also limits the range 
of accessible basin depths by the condition 

0<<f><e o . (12) 

The enumeration function in this limit corresponds to the 
logarithmic energy landscape of Debenedetti, Stillinger, 
and Shells 

uj TS ((t)) = s u-u\n(u) - (1 -u)ln(l -it), (13) 

where u — <fi/eo, 0<w<l and the superscript "TS" 
refers to the two-state model. Note that sq = is used 
in the logarithmic model in Ref. I33T 

The average basin depth is proportional to the popu- 
lation x^ s of the RS excited states in the two-state limit: 

4>* = x* TS e , x* TS = [l + e- Sa+ ^}-\ (14) 

[3 = l/k^T. The constant volume heat capacity per par- 
ticle (in fee units) is of Schottky's form22i§£ 

c^ = (/3eo) 2 a:Ts(l-a;Ts)- (15) 

This heat capacity form is continuous down to zero K 
hence, as is well known ^L^^M the two-state model 
eliminates the ideal glass transition (Tk — > 0). Thus 
variation of the Gaussian width parameters of our model 
produces the same systematic change of heat capacity 
form that Debenedetti et alm& demonstrated by ad hoc 
linear mixing of the Gaussian and binomial distribution 
CS functions. The heat capacity function for laboratory 
glasses of different fragility should then reflect the width 
of the RS Gaussians relative to the separation of their 
centers. 

The present model, which we will refer to as the two- 
Gaussian (2G) model, projects the two RS Gaussians 
onto CS enumeration function given by Eq. (|10f) . u>(<p) 
in Eq. IjlUfl is formally a linear combination of the ideal 
mixture entropy of the RS two-state model and the CS 
Gaussian term of the Gaussian model. The two terms 
are connected through the x(<p) function. x(tf>) is not a 
linear function of the two-state model [Eq. I|14|) ] when 
Cj ^ 0, although it approaches the linear limit with low- 
ering temperature (Fig. |5p.) when RS Gaussians become 
narrower (see Eq. (|2*nD . 

Both the ideal mixture term Sq(x) (first summand in 
Eq. (|10l) ) and the Gaussian term (second summand in 
Eq. IjlOl) ') are non-parabolic functions of 0. A bell-shaped 
enumeration function, which at high temperatures can be 
approximated by a Gaussian shape, is a result of com- 
bining two non-Gaussian summands in Eq. (|10f) (Fig. 
|2t>). The two-state limit of the 2G model results in a 
very asymmetric enumeration function when the excita- 
tion entropy sq is non-zero due to the cutoff of the range 
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of accessible basin energies (Eq. (|12|) : dash-dotted line in 

Pig- lib)- 

The connection between the microcanonical entropy 
to(4>) and the canonical ensemble, which permits calcu- 
lations at given temperatures, can be obtained by use of 
two thermodynamic relations^ 



and 



du>(<t>) 



N,V 



N,V 







cv ' 



(16) 



(17) 



0.2 - 



CO 




0.2 0.4 0.6 0. 



k B T/e 



FIG. 4: A (T) [Eq. JUJ] and 4>* [Eqs. <[Ei and J2TJ] vs tem- 
perature; Ai/te = 100 K, A 2 /fc B = 200 K, e /fc B = 800 K, 
so = 2.0. 

Equation (|16[) leads to the average energy of inherent 
structures 

0* = x*e - 2A , (18) 

where x* = x((p*) and 

A = (1-x*) 2 A 1 + (.t*) 2 A 2 . (19) 

The coupling parameters Ai in Eq. (|19[) are defined in 
terms of the RS distribution widths as 



of = 2fc B TA,; 



(20) 



Equations J5| and (|19|) immediately lead to Eq. (|14ll in 
the limit of narrow RS Gaussians, A; — > 0. Once Eq. 
(fT%)) is substituted into Eq. Q, one arrives at a single, 
self-consistent equation for x* which is used to obtain the 
average energy 0* in Eq. (|TH|) : 



(21) 



Here, e(x*) is the population-dependent average excita- 
tion energy 



Mechanical stability of the ideal glass state requires 

e(x*) > 0. (23) 

Equation l|21|l indicates that the introduction of fi- 
nite widths to the two RS delta functions of the two- 
state model results in self-consistency in determining the 
excited-state population which is governed by the aver- 
age ground-to-excited energy gap e(x) [cf. Eqs. lfT4*|) and 
(|2ip]. The function e(x) has a simple physical meaning. 
Going from two states with the gap e to a disordered 
material with Gaussian distributions of the RS ground 
and excited state energies makes states with lower ran- 
dom energies thermodynamically more probable. This 
effectively lowers the energy of each RS state by "sol- 
vation" energy 2Ai (Stokes shift in spectroscopic appli- 
cations). The lowering of the energy of each state is, 
however, scaled with the corresponding population and 
one gets 2xA 2 and 2(1 — x)Xi for lowering the excited and 
ground states, respectively. The dependence on ground 
and excited state populations makes the excitation en- 
ergy decrease with increasing temperature from eo + 2Ai 
at low temperature to eo — 2A 2 at high temperature (Fig. 
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e(x*) = e - 2.t*A 2 + 2(1 - x*)Xi 



(22) 



0.1 0.2 0.3 0.4 0.5 

k B T/e 

FIG. 5: cv (a) and s c (b) (in fc B units) vs temperature at 
Ai = (solid line), Ai/fc B = 50 K (dashed line), Ai/fc B = 100 
K (dash-dotted line), and Ai/fc B = 200 K (dash-double-dotted 
line); A 2 /fc B = 200 K, e /fc B = 800 K, s = 2.0, Ai and A 2 are 
assumed to be temperature independent. Numbers on plot 
(a) indicate values of Ai/fe B . 

The second thermodynamic relation [Eq. i|17fl ] allows 
one to connect the width parameter a in the CS Gaussian 
term in Eq. (|1U|) and the constant-volume heat capacity 
per molecule cy(T). The width in the Gaussian term 
can be separated into the k^T factor and an energetic 
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coupling parameter related to the RS coupling constants 
\i [Eqs. O and Hijl]: 



a 2 = 2k B TX (T), 



(24) 



The form of the CS Gaussian width in Eq. 124(1 is dic- 
tated by the classical limit of the fluctuation-dissipation 
theorem (FDT) — The energy parameter Ao(T) then 
plays the role of the trapping energy in theories of 
random-media conductivity— or the solvent reorganiza- 
tion energy in theories of electron transfer— Alterna- 
tively, Ao(T) enters the constant volume heat capacity 
obtained from Eqs. ijTU)) and (|T7|) 



c v {T) = 2/3A (T) [1 - a*{e Q - Ax* A 2 + 4(1 - x*)\i) 
where 

dx((ft) 



a 



(25) 
(26) 



The separation of the width a 2 in the k^T factor 
and Ao(T) is convenient when the latter is only weakly 
dependent on temperature. This implies that cy(T) 
is approximately a hyperbolic function of temperature 
Cv{T) oc 1/T as is empirically documented, at least for 
the constant pressure heat capacity— di* The tempera- 
ture dependence of Ao(T) [Eq. ifT^ll ] is determined by 
the gap AA = A2 — Ai in the coupling parameters be- 
tween the ground and excited states since Ao(T) has a 
sigmoidal form (Fig. 0J) with the change AA from low to 
high temperatures (Ai are assumed to be temperature in- 
dependent in Fig. A qualitatively similar sigmoidal 
form often observed in simulationsiSi 7 ^* 7 ^ 7 ^ is seen for 
the basin depth (ft*. The basin depth increases by the 
amount eo — 2AA with increasing temperature (Fig. 

The description of the glass thermodynamics in terms 
of a distribution of basin energies involves the projection 
of the whole configuration space onto a single collective 
coordinate. This projection onto the lower dimension 74 
introduces a reduced description of the system thermody- 
namics reflected by the temperature dependence of the 
moments of the collective coordinate, as is well known 
from e.g. electron transfer theory^ In fact, thermody- 
namics [Eq. (|T7Jl ] and the classical limit of the FDT [Eq. 
(|2~4*|) ] both predict that, for fluctuations caused by clas- 
sical motions, a 2 should decompose into the k^T factor 



and a coupling parameter Ao (T) . This is a significant de- 
parture from the simple temperature independent behav- 
ior of a 2 supposed to date. Such a temperature depen- 
dence modifies predictions of even the simple Gaussian 
model: no ideal glass transition occurs for temperatures 
which leave the expression sq(x*) — /3Ao(T) positive. We 
also note that the combination of disorder with the FDT 
leads to the excitation energy e(x*) [Eq. (|22|l ] decreasing 
with temperature. 

The distribution of basin energies is affected by tem- 
perature through the explicit k-gT factor in Eq. I|24|) and 
through a more complex temperature variation of Ao (T) . 
Figure EJ; shows uj((ft) at different temperatures obtained 
under the assumption that A; defined by Eq. lj2U|l are 
temperature independent. The distribution, which is al- 
most Gaussian at high temperatures, gets skewed from 
the high-energy wing at lower temperatures (cf. solid and 
dash-dotted lines in Fig.^). The low-energy wing of the 
enumeration function is not strongly affected by tempera- 
ture, and low-energy wings at different temperatures can 
approximately be brought to one master curve by a ver- 
tical shift. The high-energy wings differ, however, sub- 
stantially as temperature changes. Exactly this behavior 
of the enumeration function curves at different tempera- 
tures is reported for the 80-20 LJ binary mixture in Fig. 
4 of Ref . © 

The enumeration function taken at the average basin 
depth (ft* gives the configurational entropy (in fee units) 



s c = s (x*) - f3\ . 



(27) 



This equation predicts the existence of ideal glass tran- 
sition, s c (Tk) — 0, at a finite temperature when Ai 7^ 0. 
The Kauzmann temperature Tk tends to zero when 
Ai — > and s c > for any < x* < 1 at X x = 
(solid line in Fig. [SJd) . The constant volume heat capac- 
ity can directly be calculated from Eq. I|25|) . Calculations 
of cv(T) at constant A2 and varying Ai, both tempera- 
ture independent, are shown in Fig. |SJi. At Ai ^ the 
heat capacity drops to zero at the point of ideal glass 
transition at T — Tk > 0. There is no ideal glass transi- 
tion when Ai — 0, and the heat capacity passes through 
a broad maximum. 



III. LIQUID-LIQUID(GLASS) FIRST ORDER e(x*) remains positive at all excited state populations 
PHASE TRANSITION [Eq. Increasing the entropy of excitation allows 

one to reach the condition of vanishing Gibbs energy of 
excitation 

Depending on the value of the excitation entropy Sq 
Eq. I|21(l may have one, two, or three solutions. Only one 
solution exists for low So since the condition of mechani- 
cal stability requires that the effective excitation energy 9i x *) — e ( x *) ~ ^bTsq (28) 
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TABLE I: Best-fit parameters for model fluids (rows 1-2) and real liquids (rows 3-6). For real liquids, the fitting parameters 
are Ai, A2, to (K) and so are obtained by global fits of Eqs. 122H and 1321 to experimental heat capacities and entropies. 



Substance 


z 


T 3 


T fus 


As fus a 


eo/&B 


Ai/fc B 


A2/fcB 


so 




T K C 


80-20 BLJM d 


1 








69.3 


8.9 


26.5 


0.14 


34.35 c 


30.4 


50-50 BSSM / 


1 








78.5 


2.0 


22.3 


1.59 




14.2 


Glycerol 


8 


190 


291 


7.55 


631 


18 


32 


1.27 


136.7 


130 


Selenium 


1 


304 


494.33 


1.50 


905 


14 


27 


1.7 


210.7 


153 


Toluene 


2 


117 


178.15 


4.48 


1045 


36 


518 


5.1 


99.9 


no 9 


o-terphenyl 


2 


246 


329.4 


6.28 


2267 


45 


1133 


5.5 


204.1 


214 9 



a In fee units. 

'Kauzmann temperature (K) from extrapolation of experimental 
configurational entropy s c . 

c Kauzmann temperature (K) from the condition s c (Tk) = ob- 
tained from the 2G model. 

^Binary LJ A-B mixture with caaAb = 119-8 K, eabAaa = 
1.5, cbb/caa = 0.5, ubb/o-aa = 0.88, and density p/po = 1.2, 
po = 2.5310 28 m 3 . 

c Calculated from s c (Tk ) = in Ref. l4St Tk obtained from 
the potential energy landscape method is 34.7 K, temperature 
of vanishing diffusivity from the Vogel-Fulcher-Tammann plot is 
T = 35.97 K. 

^Binary soft sphere mixture. The reduced energy parameters from 
the fit are multiplied with caa/^b = 119.8 K for consistency with 
the 80-20 BLJM. 

9 Determined as the temperature of the first-order phase transition 
at which the entropy discontinuously drops to zero. 



in the range < x* < 1. At the temperature 

T LL = (e - AA)/s (29) 

the excitation Gibbs energy vanishes at x* = 1/2, 
(7(1/2) = 0. This is the temperature of the equilib- 
rium liquid-liquid, first-order phase transition between 
the low-temperature liquid with low concentration of ex- 
citations and the high-temperature liquid with high con- 
centration of excitationsiS*2& This transition becomes a 
liquid-glass transition when the low-temperature phase 
has its viscosity below the glass transition limit or when 
the transition occurs directly to the ideal glass state. 

The first order transition is possible for temperatures 
below the critical temperature T c and excitation en- 
tropies above the critical value sq c : 

T < T c , s > s 0c . (30) 

The critical parameters are 

k B T c = A = (Ai + A a )/2, 

e - AA (31) 

The critical temperature T c increases and the critical en- 
tropy soc decreases with enhanced disorder of the excited 
state (Fig. 0). When T < T c and s > s 0c , the equilib- 
rium transition temperature, Tll, is flanked by lower, TJ, 
and upper, T u , spinodal temperatures at which only two 
solutions of Eq. (J^IJ are possible: T ; < T LL < T u (Fig. 
EJ). 

Figure0illustrates the change in the temperature vari- 
ation of the thermodynamic parameters with increasing 
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FIG. 6: Equilibrium first-order transition temperature (solid 
lines) and the lower, T;, and upper, T u , spinodal temperatures 
(dashed lines) vs the excitation entropy so. The calculations 
are carried out for eo/^B = 800 K, Ai = 50 K, and values 
of A2/&B indicated on the plot; (soc,Xc) indicates the critical 
point. 



the excitation entropy sq. At so < %, the entropy and 
basin energy are both continuous functions of tempera- 
ture. The heat capacity passes through a broad maxi- 
mum characteristic of the two-state model and drops to 
zero at a Kauzmann temperature Tk > when Ai 7^ 0. 
At the critical excitation entropy sq — sq c , the entropy 
and energy both pass through an inflection point reflected 
by a lambda singularity in the constant volume heat ca- 
pacity. Finally, the liquid-liquid (glass) phase transition 
occurs above s 0c . The entropy drop at the transition tem- 
perature increases with increasing s to the point where 
the entropy drops to zero. This transition, as well as all 
transitions with higher entropy sq, occur directly from 
the supercooled liquid to the ideal glass state. 
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FIG. 7: Basin energy (a), configuration entropy (b), and 
constant- volume heat capacity (c) at the values of so indicated 
on the plot; Ai/Jfes = 50 K, X 2 /k B = 300 K, e /k B = 800 K. 
The values so = 3.14 and k B T/eo — 0.22 correspond to the 
critical point at which lambda singularity for the heat ca- 
pacity is seen. The dashed lines indicate metastable states 
between the lower and upper spinodal temperatures. 



IV. COMPARISON TO SIMULATIONS 

The temperature variation of the energy landscape is 
critically affected by the approximately linear temper- 
ature dependence of the width parameter a [Eq. (|24[l ] 
which contributes to the overall temperature variation of 
the distribution of basin energies 



P{4>) (x e^-^. 



(32) 



This generally non-Gaussian distribution is often approx- 
imated by a Gaussian function: 



P(4>) oc exp I — 



2r 2 



(33) 



with the empirical Gaussian width T. Computer simu- 
lations of P(<f>) and the constant volume heat capacity 
cy at varying temperature may provide insights into the 
temperature dependence of a. Approximately Gaussian 
distribution of basin energies has been found in simu- 
lations of binary LJ^^^S^ and hard sphere fiuids4^ 
Although an increase of the width with temperature is 
often seen in simulations numerical data are very 
limited. 

Temperature-dependent width T(T) can be found in 
simulations of the well-known 80-20 LJ mixtur e) 79 ! 80 by 




FIG. 8: Variance of the distribution of metabasins in the A- 
B binary LJ fluid determined from simulations by Denny et 
ffl j77 i 78 ^ an£ j ^ Doliwa and Heuej— (b). Units of energy 
and temperature are defined by the LJ energy of the A-A 
interaction potential. 



Biichner and Heuer— Also, recent extensive simula- 
tions by Denny et alSt and by Doliwa and Heuer— give 
the variance of energies of metabasins (basins separated 
by small minima which do not require activated hops at 
a given temperatur e) 72 ! 81 ) for the same system at differ- 
ent temperatures. The width of metabasin distribution 
T(T) extracted^ 8 , from the simulations by Denny et alH 
is approximately linear in T (Fig. [S^). An even steeper 
T(T) is found in simulations by Doliwa and Heuer— of 
the same system of a smaller size (Fig. |SJ)). 

The 2G model is based on four parameters: tQ, \\, 
\-2, and so- The model is tested on its ability to repro- 
duce several thermodynamic observables for a single set 
of parameters. We first apply the 2G model to simula- 
tions of model fluids and then, in Sec.0 apply it to real 
liquids. For comparison to computer experiment we fit 
the average basin energy from 2G model to combined er- 
godic parts of two cooling runs reported in Ref. for 
the 80-20 binary LJ mixture (Fig. [3*, Tabled). The fit- 
ting parameters obtained for the basin depth are used 
to calculate the configurational entropy [Eq. (|27|l ] which 
goes to zero (Fig. [3:) at the Kauzmann temperature very 
close to that obtained from simulations of s c itself and 
to To from diffusivity extrapolated to zero through the 
Vogel-Fulcher-Tammann plot (cf. columns 10 and 11 in 
Tabled). 

Since the distribution P(<f) is generally non-Gaussian, 
the empirical width T was obtained from the half- 
intensity width of P ((/))■ r 2 (solid line in Fig. EJd) is larger 
than a 2 (dashed line in Fig. |5J) since the former reflects 
the overall width arising from the ideal mixture entropy 
and the CS Gaussian term in Eq. i|10fl combined. The 
empirical width T 2 (T) rises with temperature in accord 
with the prediction of the FDT and the results of simula- 
tions. The magnitude of T 2 for the distribution of basins 
is significantly larger than the one for the distribution of 
metabasins (cf. Fig. [8] and Fig. |3Jd). However, the vari- 
ance of basin energies from MD simulations by Biichner 
and HeueriS^J. (marked "MD" in Fig. Et>) is m reason- 
able agreement with the 2G model. Finally, the constant 
volume heat capacity shows a steep rise on approach to 
the ideal glass transition, dropping to zero at Tk- This 
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FIG. 9: Inherent structure energy (a), the distribution width 
(b), configurational entropy (c), and constant volume heat 
capacity (d) for the 80-20 LJ binary mixture i 48 ' 71 ' 80 The solid 
lines refer to calculations with the present model with the 
parameters obtained from the fit of inherent structure energies 
from Ref. |^ combining data for the cooling rate 2.7 10 -4 
(squares) and 3.3310 -6 (circles). The closed points in (a) 
indicate the points which were excluded from the fit since 
they represent the loss of system's ergodicity due to finite 
cooling rate. In (b), the solid line refers to T 2 , the dashed 
line refers to a 2 , and "MD" refers to simulation results for 
r 2 by Biichner and Heuer4"*— In (d), the dashed line refers 
to the constant P heat capacity calculated under assumption 
that Xi are temperature independent. The bold solid line 
marked "MD" is obtained by numerical differentiation of the 
configurational entropy from MD simulations by Sciortino et 
al& and "2G" marks the present model. Units of energy 
and temperature are defined by the LJ energy of the A-A 
interaction potential. The parameters of the fit are listed in 
Table □ 



form of the heat capacity is supported by simulations of 
Sciortino et alml as shown in Fig. We note that the 
80-20 system was originally parameterized to represent 
the metallic Ni-P alloy 79 and that metallic glassformers 
typically show very sharp excess heat capacity functions 
relative to molecular and ionic glassformers. 83,84 

A similar fit of the 2G model to the average basin 
depth of the 50:50 soft sphere (SS) mixture reported by 
Yan et al~ is shown in Fig. ^] The parameters of the 
fit are used to calculate the distribution width T (Eq. 
iJsSJ, Fig. ITUb). the configurational entropy (Eq. (|T7|) . 
Fig. 110b). and the constant volume heat capacity (Eq. 
(|25|l . Fig. I10H). The basin energy width calculated from 
2G model (T/eAA — 0.08) is higher than the one observed 
in simulations&S (T/e AA = 0.02, marked "MD" in Fig. 

nnb). 

Despite the use of parallel tempering MD in Refs. |5l| 
and |53, the drop of the heat capacity at the peak tem- 
perature T p ~ 36.5 K (eAA/ks — 119.8 K) seen in the 
simulations is due to insufficient sampling of the phase 
space£L£2 The parameters obtained from the fit of the 
average basin depth (Table give a reasonable descrip- 
tion of the simulated cy (T) up to T p followed by a much 
stronger rise of cy(T) which drops to zero at Tk — 14 



FIG. 10: Same as in Fig. M for the 50-50 BSSM studied by 
Yan et aZi— . The parameters of the fit are listed in Table|I] In 
(b), "MD" marks the Gaussian width from parallel tempering 
MD simulations reported in Refs. IHH and l52t In (d), "2G" 
marks the present two-Gaussian model, the bold line marked 
"MC" refers to the results of Monte Carlo simulations from 
Ref. |49j. The dashed line shows the fit of the cv simulation 
data— extrapolated below the temperature of heat capacity 
drop. The functional form of cv(T) (dashed line) is from Ref. 
|H c v (T) = Ai_T Bl +A 2 T B2 , At = 2.845, Bi = -0.209, 
A 2 = 3.47710" 4 , B 2 = -5.804. 



K. This latter temperature is close to the point of van- 
ishing configurational heat capacity in the simulations. 
Note that a peak of cv(T) higher that the one reported 
m 

Refs. [H EJ andil was obtained in MC simulations 
of analogous binary SS mixture by Grigera and Parish 
within a simulation protocol outperforming parallel tem- 
pering. This implies that, once sampling is improved, 
cy(T) continues to grow beyond the drop at T p . Also 
note that extrapolation of cy(T) from the fit of simula- 
tion data by Yu and Carruzzo^ to lower temperatures 
goes even steeper (dashed line in Fig. ITUb l) than cy (T) 
from 2G model. 



V. EXPERIMENTAL CONFIGURATIONAL 
ENTROPIES AND THE KAUZMANN 
TEMPERATURE 

The configurational heat capacity of a liquid is nor- 
mally taken as the difference between the liquid and crys- 
tal entropies reported for constant pressure, although it is 
known that in many cases a part of the entropy of fusion 
is due to an increase in the vibrational entropy (arising 
from increases in the low frequency vibrational density of 
states in the liquid inherent structure a 58 i 73 i 85 ). The con- 
stant P configurational heat capacity can be calculated 
from the configurational entropy in Eq. I|27l) 



Cp = T[ df 



(34) 



The unknown parameter in this calculation is the temper- 
ature dependence of the model parameters eo, A^, and sq 
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at constant P. Spectroscopic measurements at constant 
P give "solvation energies" Ai through spectral Stokes 
shifts 8 - which are weakly temperature dependent i21 The 
Stokes shift relates to the coupling of a localized state 
to a thermal Gaussian bath. Since the defect excita- 
tions considered here may be more or less delocalized, 
it is currently unclear if the assumption (dXi / dT) p = 
is warranted. An alternative scenario might include no 
temperature dependence of the ideal glass distribution 
<7i corresponding to quenched disorder (i.e., Ai oc 1/T) 
and a standard dependence on temperature of (T 2 (i.e., 
A2 = Const). It turns out that, when the 2G model is 
applied to fit the experimental excess heat capacities of 
the liquid over the crystal, Acp = cpji q — cp. cryst , the 
results are fairly insensitive to the assumptions made re- 
garding the temperature dependence of Ai once the con- 
dition A2(T) = Const is adopted. The fit of experimental 
results is thus done with temperature-independent eo > A, , 
and so- 
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FIG. 11: (a): Configurational entropy s c vs temperature for 
liquids listed in Table U] Thick solid lines are experimental 
results and the thin solid lines are fits to the 2G model. Ex- 
perimental data stop at the glass transition temperature T g . 
Dotted lines show the extrapolation of experimental entropies 
to the zero entropy line. The fitting parameters are listed 
in Table [I] (b): Temperature dependence of the entropy of 
toluene (scaled up). The dashed lines indicate the entropies of 
metastable states terminated at the lower and upper spinodal 
temperatures. 

The fitting procedure involves simultaneous fit of Eqs. 
(|27[) and (|34|) with four fitting parameters, (Ai, A2, eo, 



and so) to experimental heat capacities Acp and exper- 
imental configurational entropies >^ The range of energy 
parameters is restricted by the condition of mechanical 
stability of the ideal glass state [Eq. jjZEt ] ■ The configu- 
rational entropy at constant P can be determined exper- 
imentally from the entropy of fusion Asf us and Acp(T) 
(both in fee units) 



s c (T) = As {us + I (Ac P (T')/T')dT' 



(35) 



The 2G model outlined in Sec. [H] assumes that each 
molecule represents one excitable unit. While this is true 
for atomic glasses like selenium, for more complex com- 
pounds one needs to introduce the number z of indepen- 
dently excitable (i.e. rearrangeable) states per molecule 
or formula unit^S The parameter z is taken from Takeda 
et al£& (TableP) and is used to multiply the heat capacity 
in Eq. Q-i'lH in fitting the experimental data. The results 
of the fit for four glassformers are listed in Table [I] 
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FIG. 12: Comparison of the experimental values of Acp (solid 
lines) with Acp calculated from the 2G model (dashed lines) 
by simultaneous fit to the experimental s c and Acp data 
above T g . 

The examination of Table [I] shows that Ai > for 
all fluid studied, indicating that Tk > 0. The distri- 
bution of excited states is much narrower in the case 
of covalent and hydrogen-bonded liquids (selenium and 
glycerol) compared to molecular liquids (toluene and o- 
tcrphcnyl). In contrast, the RS distribution of the ideal 
glass is almost invariant among different glassformers. 
When A 2 » Ai the configurational entropy s c (T) gains 
a bend close to the Kauzmann point (Fig. Ill|) result- 
ing in the actual Tk from s c (Tk) = smaller than the 
corresponding value from extrapolation of experimental 
entropies (e.g., selenium in Table P). 

The most interesting result of our analysis is the 
low-temperature behavior of fragile molecular glasses 
(toluene and o-terphenyl) . These substances are char- 
acterized by high disorder of the excited state (A2 3> Ai, 
Fig. fH3|) and, in addition, high entropy of excitation (Ta- 
ble P). It also turns out that sq from the fit is higher 
than the critical excitation entropy which is close to 2.0 
for both liquids. The fact that sq is more than twice 



12 



higher than sq c ensures low first-order transition tem- 
perature, well below the critical temperature (543 K for 
toluene and 1156 K for o-terphenyl) . The first-order tran- 
sition temperature in fact falls in the unobservable range 
between Tk and T g where the entropy discontinuously 
drops to zero producing a similar drop in the heat ca- 
pacity fFies. fTTI and I12p. It may be therefore suggested 
that fragile liquids resolve the Kauzmann paradox by a 
first order liquid-glass transition. We note that both for 
toluene and o-tcrphenyl the lower spinodal temperature 
is below the point when metastable entropy crosses the 
zero entropy line while the upper spinodal temperature 
is above T g for toluene and almost coincides with T g for 
o-terphenyl. Since the first order transition is below T g , 
the equilibrium passage along the solid line in Fig. ITTb 
is unlikely thus suggesting hysteresis of the heat capacity 
between the cooling and heating runs. 




e/k B ,K 

FIG. 13: Real-space distribution of energies for selenium 
(solid lines), glycerol (dashed lines), and toluene (dash-dotted 
lines) at their corresponding Kauzmann temperatures. The 
parameters eo and Xi used in the calculations are taken from 
Table |H 

The value So, which remains an empirical parameter 
of the 2G model, can be compared to the entropy cost 
of creating a density wave in density-functional theories 
of aperiodic struct ures ffiSS ' The average entropy of the 
"entropy droplet" in Wolynes's mosaic model is 

3 5 
s = - In (arl/ir) - -, (36) 

where a represents the rms displacement from the lattice 
site and ro is the mean lattice spacing. Invoking the 
Lindemann ratfoiS a 1 ^ 2 ^ = 10, one gets sq = 2.7, which 
falls in between entropies for strong and fragile liquids in 
Table D 

Caution is needed in these interpretations since the 
model is in the early stages of evaluation and there 
are four parameters even for simply constituted glasses 
(z = 1). One of these parameters may be disposable. It 
is apparent from Table |U and Fig. ^3 that the ground 
RS Gaussian width best fitting the various data, while 
non-zero (as expected for a non-crystalline ground state), 
is small relative to the excited state Gaussian (except 



for the stronger liquids) and not varying much between 
the different systems. It could probably be given a fixed 
value, reducing the disposable parameters to 3 for simple 
glasses and 4 for flexible molecule glasses, where the 4th 
parameter can be fixed from molecular considerations; 88 



VI. CONCLUDING REMARKS 

We have shown that by introducing a realistic form 
for defect-like excitations in glasses, the basic "excita- 
tions" model of the glass transition can be developed in a 
form that bridges the gap between previous over-simple 
models and the random energy model of Derrida4°» In 
other words, we have provided a physical basis for the 
previously empirical "logarithmically modified Gaussian" 
model of Debenedetti et alm& The model is fundamen- 
tally non-Gaussian in configuration space. It recognizes 
the role of fluctuations within the FDT in making the 
landscape temperature-dependent as evidenced by the 
behavior of (meta)basin energy variances from molecu- 
lar dynamics simulations ^M^L^S^SL 

The model predicts a possibility of first-order liquid- 
liquid (glass) transition when the entropy of excitations 
exceeds its critical value and the temperature falls below 
the critical point. For fragile liquids characterized by a 
broad distribution of excitation energies and high entropy 
change per excitation the transition temperature is low. 
While most known liquid-liquid transitions for strong liq- 
uids are at high temperatures^ the observation of such 
a transition for fragile supercooled triphenyl phosphite 76 
supports this trend. The fit of the model to experimen- 
tal entropies and heat capacities of fragile toluene and o- 
terphenyl results in the first-order liquid/ideal glass tran- 
sition between T g and the experimental Kauzmann tem- 
perature. It seems therefore reasonable to suggest that 
fragile liquids release the excess entropy by a first-order 
transition to the glassy state. 

The present model belongs to a class of mean-field two- 
state models in which the average excitation energy drops 
linearly with the increase in the population of the ex- 
cited state [Eqs. (12111 and 1)22(1 ], Negative excitation en- 
ergies are prohibited by the condition of mechanical sta- 
bility, and crossing the zero point of the excitation Gibbs 
energy is driven by the excitation entropy the magni- 
tude of which is correlated with glass fragility. 90 Another 
physical realization of this model is the coupling between 
molecular excited states through long-range interactions. 
In case of optical excitations of molecules coupled by 
long-range dipolar forces the change in the excitation en- 
ergy is realized through the reaction field proportional 
to the number of excited molecules. A mean-field de- 
scription, mathematically equivalent to the present 2G 
model, then results in transition to excitonic condensate 
in molecules coupled through their transition dipolesSi 
or to a non-polar/paraelectric phase transition in dipolar 
two-state fluids^ In the present model, disorder is re- 
sponsible for the trapping energy playing the role of the 
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